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ABSTRACT 

We present high resolution 8 km s~^ ) spectra of the QSO Q0201+365 
obtained with HIRES, the echelle spectrograph on the 10m W.M. Keck 
Telescope. Although we identify over 80% of the absorption features and 
analyze several of the more complex metal-line systems, we focus our analysis 
on the damped Lya system at z = 2.462. Ionization simulations suggest the 
hydrogen in this system is significantly neutral and all of the observed metals 
are predominantly singly ionized. We measure accurate abundances for Fe, Cr, 
Si, Ni and place a lower limit on the abundance of Zn: [Fe/H] = —0.830 ± 0.051, 
[Cr/H] = -0.902 ± 0.064, [Si/H] = -0.376 ± 0.052, [Ni/H] = -1.002 ±0.054 
and [Zn/H] > —0.562 ± 0.064. We give evidence suggesting the actual Zn 
abundance is [Zn/H] ^ —0.262, implying the highest metallicity observed at 
a redshift z > 2. The relative abundances of these elements remains constant 
over essentially the entire system (f» 150 km s~^ in velocity space), suggesting 
it is well mixed. Furthermore, we use the lack of abundance variations to infer 
properties of the dust responsible for element depletion. Finally, we discuss the 
kinematic characteristics of this damped Lya system, comparing and contrasting 
it with other systems. The low-ion line profiles span ^ 200 km s~^ in velocity 
space and have an asymmetric shape with the strongest feature on the red edge. 
These kinematic characteristics are consistent with a rotating disk model. 



Subject headings: galaxies: evolution — galaxies: formation — quasars: quasars — 
absorption lines — quasars: individual (Q0201+365) — galaxies: abundances 
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1. INTRODUCTION 

This paper is the second in a series devoted to studying the metal content of 
high-redshift galaxies and their progenitors. Our primary objectives are 

(1) to record the emergence of metals in galaxies, 

(2) to trace the mean cosmic metallicity from z ~ 4.5 to the present, 

(3) to determine the kinematic state of galaxies from z ^ 4.5 to the present. 

We are implementing this study using HIRES, the echelle spectrograph on the Keck 10m 
telescope ( |Vogt 199^) , to obtain high- resolution spectra of QSOs with foreground damped 
Lya systems. The damped Lya systems are a population of neutral gas layers exhibiting 
properties indicating they are either galaxy progenitors, or well formed galaxies detected 
during an early evolutionary phase. Recent studies indicate the comoving density of neutral 
gas in damped systems at 2; ~ 3.3 is comparable to the density of visible stars in current 
galaxies. At lower redshifts, the comoving density of neutral gas decreases with time in a 
manner consistent with gas consumption by star formation ( [Wolfe et al. 1995|) . Therefore, 



studies of the metal content of the damped Lya systems enable one to trace the chemical 
evolution of representative galaxies from a presumably metal-poor gaseous progenitor phase 
to metal rich epochs when most of the baryons are in stars. As a result, the age-metallicity 
relation, kinematic conditions, etc., deduced from the damped Lya systems should tell us 
more about the history of galaxies at large redshifts than analogous relations deduced from 
old stars found in the solar neighborhood (jEvardsson et al. 19931). 



In a previous paper we presented echelle spectra of the z = 2.309 damped system 
toward PHL 957 at a spectral resolution of ~ 8 km s~^ (FWHM) and signal-to- noise ratio 
of ~ 35:1 ( [Wolfe et al. 19941) . By fitting multiple Voigt velocity components to low-ion 
transitions such as Zn II 2026, Ni II 1741, and Cr II 2062 we obtained accurate abundances 
for Zn, Ni, and Cr in the neutral gas. The Zn and Cr abundances were accurate because we 
resolved Zn II 2062.664 from Cr II 2062.234 for the first time in a QSO absorption system. 
We found that the abundances relative to solar were low: [Zn/H] = — 1.55±0.11, [Cr/H] = 
— 1.79±0.10, and [Ni/H] = — 2.13±0.08. The Zn abundance is especially significant because 
Zn is relatively undepleted by grains in the ISM of the Galaxy ([Sembach et al. 1995|) and is 



presumably unaffected by dust which may be present in damped Lya systems ([Fall and Pei 



1993[) . We also found the line profiles to be asymmetric in the sense that the low column 
density gas was found in absorption only at velocities higher than the high column- density 
gas. The kinematics can be explained by the passage of the line of sight through a rotating 
disk in which the density of clouds decreases with radius and with perpendicular distance 
from midplane. 
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The purpose of this paper is to present HIRES spectra for Q0201+365, & V = 17.5 
QSO with emission redshift Zem = 2.49. While we identify more than 80% of the absorption 
features and find a total of 13 metal-line redshift complexes, the focus of this paper is on 
the damped Lya system at z = 2.462 ( Lanzetta et al. 1991 , Lu et al. 1993| ). Lu et al. 1993 



studied this system at a resolution of ~ 50 km s . They fitted a Voigt damping profile to 
the Lya absorption trough and found A^(HI) = 2.4 x 10^°cm~^. They also identified the 
metal transitions Si II 1190, Si II 1193, Si III 1206, Si II 1260, and Fe II 1144. Because 
these transitions are (a) saturated and (b) in the Lya forest, neither the abundances nor 
kinematics were accurately measured. The present study represents a major improvement 
over the previous work for the following reasons. First, we obtain spectra at a resolution 
of ~ 8 km s~^ and a typical signal-to-noise ratio of 33:1. Second, in contrast to the 
previous work, we focus on metal lines redward of Lya emission where confusion with Lya 
forest lines is absent. Third, because of the higher accuracy of the data we focus on weak 
unsaturated transitions of ions expected to dominate the ionization state of gas in neutral 
clouds. In fact we establish accurate element abundances for Fe, Si, Ni, and Cr, and a lower 
limit on the abundance of Zn. Furthermore, we use computer simulations to investigate the 
ionization of this system. We also analyze the relative metal abundances and comment on 
the characteristics of dust grains in this system. In addition, we examine a Lya absorption 
system at z = 1.955 and use kinematic and abundance arguments to suggest it may be a 
damped Lya system. Finally, we discuss the kinematics of the 2;=2.462 damped absorption 
system, contrasting its features with several other systems toward Q020 1-1-365, as well as 
other damped Lya systems measured with HIRES ( [Wolfe et al. 199^ ). 

The paper is organized as follows. In § 2 we describe the data acquisition and reduction 
techniques, present the spectra (Figure |l|) and give a nearly complete absorption line list in 
Table 0. We detail the analytic methods utilized throughout the paper in § 3. In § 4 we 
present velocity profiles of the most significant metal line systems along with the VPFIT 
package solutions where applicable. § 5 presents the ionic column densities of the two 
systems associated with the damped Lya profile at ^=2.46. In § 6 we argue that the degree 
of photoionization of the damped Lya system is low. § 7 gives the results of the abundance 
measurements and discusses the possible depletion of the gas-phase metals by dust grains. 
We also describe the kinematics of the Lya system. Finally, § 8 summarizes the results and 
gives concluding remarks. 
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2. DATA 

In this section we present the HIRES spectra, detailing the techniques used for the 
acquisition and reduction of the data. Table |^ gives an absorption line list with measured 
equivalent widths and la errors and identifies over 80% (> 88% redward of the Lja forest) 
of the features. 



2.1. Acquisition 

We observed Q0201+365 with the HIRES echelle spectrograph on the 10m W.M. Keck 
Telescope on three separate nights for a total integration time of 9.7 hrs. Table |l] presents 
a journal of the observation dates, exposure times, wavelength coverages and resolution 
of the data. Unfortunately, the signal to noise ratio (SNR) of the data was significantly 
limited by clouds. We used the kv38 filter to block out 2nd order blue wavelengths, the C5 
decker plate with 1.1" slit, and standard 2x1 binning on the 2048 x 2048 Tektronix CCD. 
This setup afforded a resolution ranging from 7.2 — 8.0 km s~^ and wavelength coverage 
from 4720 — 7180 A over 26 orders. Gaps are evident between orders redward of ~ 5100 A 
where the free spectral range of the echelle exceeds the format of the CCD. For reduction 
and calibration we took 300s exposures of the nearby standard star BD+28411 and images 
of quartz and Th-Ar arc lamps. 

2.2. Reduction 

The 2-D CCD images were reduced to 1-D uncalibrated spectra with a software 
package kindly provided by T. Barlow (1995). In short, the package subtracts the baseline 
and bias from the 2-D frames, determines the gain from dark images converting digital 
numbers (DN) to e~ counts, and extracts the data from the 2-D images by tracing the 
bright standard-star profile. It also performs sky subtraction and the removal of cosmic ray 
events. 

We wavelength calibrated the data of each night by fitting low-order Legendre 
polynomials to our Th-Ar arc frames, properly correcting to vacuum heliocentric 
wavelengths. We obtained continuum fits to the spectra by fitting high-order Legendre 
polynomials to the spectral flux in each order. The continuum fits were then used 
to normalize the flux to unity. We also calculated a la error array assuming Poisson 
statistics, ignoring the errors associated with sky subtraction and continuum fitting (i.e. 
a = J Noh} + Asky). Finally we coadded the spectra, weighting by the calculated SNR of 
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each image while rejecting spurious values assumed to arise from cosmic rays or poor sky 
subtraction. The final average SNR is ~ 30 with the exception of Order 75 (SNR ^ 22) 
where only one night of data was acquired due to a difference in the alignment of the CCD 
on the first night. The spectra are presented along with Icr error arrays in Figure |I|. 

2.3. Absorption Lines 

Table ^ lists the wavelengths, equivalent widths and la errors for all absorption line 
features which exceed the 5a limit in equivalent width as measured by techniques similar to 
those of Lanzetta et al. (1991). We believe 5a is a conservative but appropriate limit for 
this data. The reported wavelengths represent rough estimates of the centroids of complex 
line profiles and should be taken only as guides for differentiating between features. Because 
almost all of the transitions are resolved, an accurate determination of the wavelength of 
every feature would require months of laborious Voigt profile fitting. This task bears little 
scientific merit for the task of the present paper which is to examine specific systems. For 
this reason, we carried out detailed profile fits for only a selected subset of the lines in 
Table |^. We do not report equivalent widths for those features identified as sky absorption 
lines or complicated multi-system blends. 

Table ^ also includes the transition names and approximate redshifts of those features 
we successfully identified. Identification proceeds in a largely ad hoc fashion, with the 
initial emphasis placed on finding C IV and Si IV doublets. Once we composed a list of 
redshifts for the metal-line systems, we attempted to match the remaining features with 
the strongest metal-line transitions. Finally, we compared the line profiles of the transitions 
of a given redshift system in velocity space for conclusive identification. Although fairly 
laborious, this approach is highly effective and essentially error-proof. 

By comparing the object frames with the similarly reduced standard star images, we 
identified night-sky emission and absorption features in the spectra. These features are 
labeled appropriately in Table || along with all other identified spurious features. 

Table ^ lists the rest wavelengths and oscillator strengths for all of the metal transitions 
analyzed in this paper. Almost all of the values are taken from Morton (1991). It is 
important to be very clear what values are assumed in order to allow for consistent 
abundance comparisons. 
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3. ANALYTIC METHODS 

This section describes the least squares hne profile fitting method and the apparent 
optical depth method used to analyze our metal line systems. 

3.1. VPFIT (Least-Squares Line-Profile Fitting Method) 

With the aid of the VPFIT fitting package kindly provided by R.F. Carswell we 
performed least-squares fits of Voigt profiles, produced by individual Gaussian components, 
to our metal lines for several of the absorption systems toward Q020 1+365. The VPFIT 
package fits Voigt profiles to an absorption system, simultaneously determining the redshift, 
column density, b values (where the Doppler velocity b and velocity dispersion a are related 
by 6 = V^cr) and associated errors of the individual components while minimizing a 
parameter matrix . When performing a fit, one can tie several transitions together forcing 
the redshift and b values of the different transitions to match while allowing the column 
densities to vary individually. 

In performing the fits we first isolated the broadest resolved transition in a given 
absorption complex. We found we could always obtain a reasonably accurate profile fit to 
this transition. We then tied it together with the other associated transitions. Given the 
inherent differences between low and high-ion line profiles, we chose to treat the two types 
separately. If necessary, we added or removed a velocity component to those transitions with 
line profiles that have features not evident in the other transitions. In several cases, this 
significantly lowered the final value. In all our fits we assumed bulk motion dominates 
thermal motion because damped Lya systems are relatively cool (T < 1000 K) and the 
transitions arise in metals that are comparatively massive. We also set a minimum value 
for the b parameter at 3 km s"^ to prevent the package from fitting features narrower than 
4 pixels, the FWHM of the line spread function. 

3.2. Appeirent Optical Depth Method and Hidden Component Aneilysis 

Savage and Sembach (1991) have stressed that measuring column densities with the 
line profile method does not always account for hidden saturated components. These 
saturated components may be under represented in the line profile analysis, leading to 
significant errors in the measured ionic column densities. Therefore, as a check on our 
VPFIT results we performed a hidden component analysis of those ions where multiple 
transitions were observed (e.g. Ni"*" for the damped Lya system a,t z — 2.462). The analysis 
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involves calculating Na{v), the apparent column density per unit velocity, for each pixel 
from the optical depth equation 

iv„(.) = (1) 

where Ta{v) = ln[/j(t>)//a(t>)], f is the oscillator strength, A is the rest wavelength and Jj 
and la are the incident and measured intensity. Comparing Na{v) deduced from two or 
more transitions of the same ion, one finds the stronger transition will have smaller values 
of Na{v) in those features where hidden saturation is present. Thus, one can ascertain 
the likelihood of saturated components for ions with multiple transitions. Furthermore, 
summing Na{v) over the proper velocity intervals serves as an excellent check on ionic 
column densities measured with the VPFIT package, and as a further check on the existence 
of saturation. 



4. VELOCITY PLOTS AND VPFIT SOLUTIONS 

This section presents the velocity profiles of the most complex metal line absorption 
systems toward Q0201+365. For several of the systems, we superimpose the solutions of 
our least square fits from the VPFIT package. In all plots a dashed vertical line is drawn 
for reference, usually identifying the strongest feature present. For clarity, we have plotted 
features not related to the given system with dotted lines. 



4.1. Damped Lya System (z = 2.462) 

Figure ^ presents the line profiles and fits of the transitions associated with the damped 
Lya system at 2;=2.4620. The velocity centroids for the Gaussian components of the fit to 
the low-ion profiles are denoted by the short vertical lines above the Fe II 1608 profile. The 
velocity centroids for the high-ion fits are the vertical lines above the Si IV 1393 profile. 
We found 13 individual components were required for an optimal fit to the low-ions. The 
multi-component fit to the transitions Fe II 1608, Ni II 1741, Ni II 1751, Cr II 2056 and Si 
II 1808 has a reduced xt = 1-05 with a probability P^2 = 0.216. Table Hal lists the redshift. 



b value and column density along with la error of every velocity component in the VPFIT 



^ Note there are two 'subsystems' within the velocity space of the damped Lya profile: the present 
subsystem and the z=2.457 subsystem discussed in §4.2 
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solution. The quality of the fit reflects the high degree to which the low-ion profiles track 
one another. 

On the other hand, 10 components (only 8 over the same region as the low ion profiles) 
were necessary for an optimal fit to the high-ion transitions C IV 1550 and Si IV 1393. 
These components are significantly broader (higher b values) than those of the low-ions. 
This fit is not as accurate as the low-ion fit [xl = 1.805) because it is much more difficult 
to fit broad shallow components such as those around ~ —260 km s~^ . The results of this 
fit are presented in Table 

Figure ^ shows velocity profiles of 4 transitions we chose not to fit with the VPFIT 
package. Because the Al II 1670 profile is so highly saturated, we were unable to accurately 
fit it with the other low-ion profiles. In particular, the VPFIT package could not properly 
model the heavy absorption at ^ —140 km s~^ . The Al III 1862 profile exhibits 
characteristics of both the low and high-ion transitions and therefore was difficult to fit it to 
either. Instead, we measured its column density with the apparent optical depth method. 
Finally, the 3 transitions Zn II 2026, Cr II 2062 and Zn II 2062 all suffer from significant 
blending: Zn II 2026 is blended with both Fe II 2374 and Fe II 2600 from 2 other systems, 
and Zn II 2062 and Cr II 2062 are blended with each other. As a result all three transitions 
were excluded from the fit. 



4.2. Companion System at z=2.457 

The velocity profiles and least-square fits of the other subsystem (found at z=2.457) 
associated with the damped Lya system at 2;=2.46 is shown in Figure ^ A 12 component 
fit to the low-ion transitions Fe II 1608, Si II 1526 and Al II 1670 was optimal and yielded 
a reduced xl = 1-50. Difficulties arose in fitting the Fe II 1608 profile because the second 
component appears to have a significantly lower b value than the same component in the 
other low-ion profiles. This may be a result of thermal broadening. We also fitted the 
high-ions SilV 1393, 1402 and CIV 1548 and again found broader, fewer components were 
required in the best fit. Tables and [SB list the results of the fits. 



4.3. Possible Damped System at z=1.955 

Figure ^ presents the velocity profiles for a redshift system with 2;= 1.955 which is a 
possible damped Lya system. We chose not to fit this system because significant blending 
with other systems prevented us from determining which features are associated only 
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with the 2;=1.955 system. In particular, we have been unable to determine whether the 
absorption seen at positive relative velocities in Figure ^is due to the z=1.955 system. 

Lu et al. (1994) could not confidently classify the system as damped because the noisy 
Lya profile could not be fitted with a Voigt damping profile. The metal lines have many 
resolved components and exhibit a large velocity interval (~ 250 km s"^ ). Figure ^ plots 
the observed Lya profile, a Voigt damped Lya profile, and the Fe II 1608 line profile. The 
Voigt Lya profile was derived by letting iV(HI) = 1.5 x 10^° cm~^ at z=1.955. It appears 
the low-ions satisfy the metal criteria associated with damped Lya systems; i.e., the low-ion 
metal profile is significantly narrower than the Lya profile and its velocity centroid is near 
that of the Lya profile (Wolfe et al. 1993). An apparent optical depth measurement of the 
Cr II and Si II ions gives iV(Cr+) = 12.73 ± .070 and A^(Si+) = 15.25 ± 0.03 suggesting a 
lower limit of log[A^(H)] > 19.0 and 19.7 respectively, assuming cosmic abundances and no 
dust depletion. This analysis indicates the z= 1.955 system may be damped. Furthermore, 
Figure |^ reveals evidence for C I 1656, a transition previously seen only in damped Lya 
systems. However, even with our high resolution data the proper classification remains 
inconclusive. 



4.4. System at ^=2.325 

The velocity profiles and VPFIT profile solutions for the absorption system at 2;=2.325 
are shown in Figure ^ Note the feature at —80 km s~^ in the CIV 1548 transition is due 
to a blend with the CIV 1550 transition from the system at 2;=2.320. By contrast with the 
damped system, the Al III 1854 transition was best fitted with the low-ion solution. We 
tied the Si II 1526, Fe II 1608, Al II 1670, and Al III 1854 transitions together and fitted 
the CIV doublet separately. Tables ^ and ^ present the results of the 2 fits. Unlike the 
damped system, the high and intermediate ions more closely follow the low-ion profiles. 
Because the high-ion profiles have several components with higher b values, however, it is 
still impossible to fit them together. 

Figure || plots the Si II 1526, C IV 1550 and Lya profiles for the z=2.325 system. All 
three transitions span nearly the same velocity interval and track one another moderately 
well. As expected the Lya profile tracks the low-ion profiles (in particular the component 
at f ~ 20 km s~^ ) more closely than the high- ion profiles. Although we did not fit the Lya 
profile. Figure ^ is consistent with the profiles of Lyman limit systems. 
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4.5. Mg II Systems at 2=1.476, 1.699 

Figure ^ shows the velocity profiles of two Fe II and Mg II transitions for metal line 
systems found at the redshifts z=lA76 (a) and ^=1.699 (b). Both sets of transitions are 
relatively complex and most of the absorption spans moderate velocity intervals (~ 100 
km s~^ ). However, the Fe II 2374 and Fe II 2586 transitions in the z = 1.476 system may 
exhibit an additional feature at ~ —275 km s~^ . Because the predicted wavelengths for the 
other Fe II lines fall in the inter order gaps, we cannot verify the reality of this feature. 

Figures 10a and 10b are hidden component analyses of the Fe transitions from the two 
systems. Fe II 2586 was not plotted in either hidden component analysis because it closely 
traces the Fe II 2374 profile in each system and would clutter the figures. In Figure 10a, the 
weakest Fe transition (Fe II 2374) has significantly larger Na{v) values, while the strongest 
Fe transition (Fe II 2600) has the smallest Na{v) profile. This is direct evidence of hidden 
saturated components. In fact, because of these hidden components, we found it impossible 
to fit all of the transitions together in the 2;= 1.476 system. We were successful in fitting the 
Fe II 2344, 2586 transitions, however, and the results are presented in Table 0. 

There is little evidence of hidden saturated components in the 2;=1.699 system. What 
the hidden component analysis does reveal, however, is a blend in the Fe II 2374 profile 
with Al III 1854 from the z = 2.457 at velocities greater than pa km s~^ and a blend in 
the Fe II 2600 profile with Fe II 2374 from the z = 1.955 system at velocities lower than ^ 
—70 km s^^ . Unfortunately, these blends, and the low SNR in the Fe II 2586 and Fe II 
2374 profiles prevented accurate profile fits. 

5. IONIC COLUMN DENSITIES 

This section presents the ionic column densities of the two systems associated with the 
damped Lya profile at z = 2.46. We perform hidden component analyses where applicable, 
and use both the line profile and apparent optical depth methods to measure column 
densities. 



5.1. Damped Lya System at z = 2.462 



The hidden component analysis of the Ni II 1741, 1751 transitions for the damped 
Lya system at z = 2.462 is shown in Figure |Tl|. With a few minor exceptions, the Na{v) 
curves for the two Ni II transitions match within la, suggesting no hidden saturated 
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components. Therefore, abundances based on column densities inferred from profile fitting 
or the apparent optical depth method should not suffer from significant hidden saturation 
effects. Table ^ lists the measured ionic column densities and la errors for Fe"*", Cr"*", Si"*", 
Al"*" and Ni"*" for the damped Lya system at z=2.462 as measured by both the line profile 
(VPFIT) and apparent optical depth methods. For the line profile method we summed 
the column densities of the individual components of each transition and calculated the la 
error in the total value with standard least squares techniques. 

We adopt a final value for the ionic column density for all of our measurements 
by averaging the two values and adopting the VPFIT errors. We found the apparent 
optical depth method underestimates the error, particularly transitions which are nearly 
saturated (e.g. Fell 1608). For nearly saturated or very weak lines, errors associated with 
sky subtraction and continuum fitting will play a significant role, yet these errors are not 
included in the la array. The VPFIT package estimates errors based on the fit of all profiles 
and is influenced more by the deviation of the data from the fit. Therefore, we chose to 
adopt the VPFIT errors for all of the ions. Note that in almost every case the calculated 
values from the two methods match, further indicating no hidden saturated components. 

As noted above, we can place only a lower limit on the Zn"*" column density because of 
blending. We find the Zn II 2026 transition to be dominated by Fe II 2374 associated with 
the z = 1.955 absorption system, and Zn II 2062 to be partially blended with the stronger 
Cr II 2062 transition. Although Zn II 2026 cannot be extracted from the Fe II profile, the 
Zn II 2062 profile can be used to estimate a lower limit on the Zn"*" abundance. Figure |l^ is 
a hidden component analysis of the Cr II 2056 and Cr II 2062 transitions. Note that the 
Cr II 2062 (dotted) profile dominates the Cr II 2056 profile over the entire velocity interval, 
as expected if Zn II 2062 is present. This analysis reveals two features around 40 and 60 
km s~^ evident only in the Cr II 2062 profile. We suggest that these two components are 
not due to Cr+ absorption, but are components of Zn II 2062 at the high velocity edge of 
the profile. In the velocity space of Zn II 2062 these components correspond exactly to 
the strongest components of the other low-ion profiles at z=2. 46258 and z=2. 46280 (i.e. 
at V = —20 and km s~^ in the Zn II velocity space). Our hidden component analysis of 
Cr II reveals a third component of Zn II at 7139.66 A (—60 km s~^ in Figure |12D, also 
corresponding to a significant feature in the other low-ion profiles. 

Summing the column densities over velocity space with the apparent optical depth 
method, we find we find logiV(Zn+). = 11.753 ± 0.116, 12.087 ± 0.066, and 12.060 ± 0.074 
respectively for the three components (blue to red). The resulting lower limit for the total 
Zn"*" abundance, log A^(Zn"'")^ = 12.468 ± 0.046 where the error reported is only useful as an 
indication of the uncertainty in this lower limit value. 
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Table |^ also lists ionic column densities for the s-process transitions Pb II 1682 and Ge 
II 1602. Because these transitions are so weak, we the used linear curve of growth to infer the 
column densities from the respective rest-frame equivalent widths. We report the following 
values for the rest-frame equivalent width of the transitions: H^(Pb) = 1.78 ±0.47 x 10~^ A 
and W{Ge) = —3.4 ± 4.6 x 10"^ A. The negative equivalent width value for Ge II 1602 
is certainly a function of the continuum error and is consistent with a null detection (less 
than even a la detection). We have chosen, therefore, to report its column density and 
abundance in terms of a 3a upper limit. The Pb II 1682 value is significant at the 3a level, 
but could also be explained through a sizable continuum error, improper sky subtraction, or 
an unidentified blend. Although Sn II 1400 and Ga II 1414 both lie within our wavelength 
coverage, they are overwhelmed by blends with transitions from other metal line systems 
and could not be analyzed. 

5.2. Companion Subsystem at z=2.457 

Table ^ presents the measured ionic column densities for the transitions of the 
companion system at z=2A57. Final values are adopted according to the criteria described 
above. In nearly all of the low-ion transitions the measured column densities are ~ 20 times 
lower than those from the z=2.462 subsystem. On the other hand the column densities of 
the high-ions are nearly the same, possibly indicating this system is in a higher state of 
ionization, that the two ion groups are kinematically disjoint, or that the same high-ion gas 
envelopes two dissimilar low-ion configurations. 

6. IONIZATION 

This section investigates the photoionization of the damped Lya system at z=2.462. 

6.1. Lya Profile 

Figure is a velocity plot of the Lja profile and Fe II 1608 profiles associated with 
the two absorption systems at z=2.462 and z=2.457 together with the velocity profile of the 
Fe II 1608 transitions. Unrelated features at ~ —800 and 500 km s~^ have been removed for 
clarity. We derived the Voigt profile by distributing the total HI column density, logA^(H°) 
= 20.38, into the 23 Fe II 1608 velocity components weighted by their corresponding column 
density fractions {N{Fe'^)-/N{Fe'^)^^^^). This treatment is proper provided (a) the Fe II 
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1608 line profiles accurately trace the HI gas, as one would predict for a sufficiently neutral 
system, and (b) [Fe/H] is constant across the entire velocity profile, which is predicted for a 
well mixed system. As Figure demonstrates, the resulting Voigt profile is well fitted to 
the low resolution data. Table ^ lists the column density, redshift and b values of the 23 
components as adopted. Here, about 5% of the HI gas is located in the companion system 
at z=2.457. We found one could place no more than 15% of the HI gas in the companion 
system before significantly distorting the left wing of the Lya line profile. 

6.2. Ionization Models 

6.2.1. Neutral Hydrogen Model 

Table |10| shows that all but the most abundant components of our damped system 
have log[A^(HI)] < 19.5. If these clouds were isolated structures in the IGM, they would 
be highly ionized. This would markedly reduce the accuracy of abundance determinations 
based on the assumption that most of the hydrogen is neutral and most of the metals are 
singly ionized. 

In order to address this problem we ran several ionization simulations with the aid 
of a program developed by Vincent Virgilio to investigate the predicted neutral hydrogen 
fraction (H'' / H) in the damped system. Figure 14a shows the neutral hydrogen fraction 
plotted against log[A^(H°)] measured from either face of a plane-parallel layer with constant 
H volume density, uh = 0.1 cm~^. The layer is subjected on both sides to attenuated 
power law continuum radiation as calculated by Madau (1992) for a redshift of 2=2.46 
with a mean intensity of J^, = .195 x 10~^^ ergs s^^ cm~^ Hz^^ sr^^ at 1 Rydberg. The 
simulation assumes a temperature of 10^ K (purely for determining the recombination 
rate of H"*") while satisfying the ionization and transfer equations in a large number (100) 
of parallel discrete cells. Each face of the plane parallel layer is illuminated by uniform 
radiation with incidence angles covering 27r steradians. The calculation also assumes zero 
source function (i.e. no reemission). Although this model is rather simplified, we find 
similar CLOUDY (version 84.12: Ferland 1991) calculations are in good agreement with 
our results. Figure |15| compares the neutral fraction predictions of our model with the 
corresponding predictions by CLOUDY. The 2-sided illumination simulation predicts a 
lower neutral fraction than CLOUDY even though CLOUDY assumes only perpendicular 
incidence. Thus, our simulations reveal that the 1-sided assumptions inherent in CLOUDY 
are probably underestimating the degree of ionization in optically thick absorption systems. 

Our simulations predict a uniform layer with A^(H'^) = 2.4 x 10^° cm^^ and 
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uh = 0.1 cm~^ will be 96% neutral. The results should be similar for a more realistic 
multicomponent system, since we expect the majority of the 23 components in the z=2.46 
system to be shielded from ionizing radiation by gas both above and below the location of 
the cloud in the layer, thereby maintaining neutrality (Figure 14b). Therefore, we measure 
the metal abundances of our damped Lya system by assuming the abundance of element X 
to equal A^(X"'")/A^(H°), where is the first ionization state of element X. Note, however, 
that this result is very dependent on the value oi 71^. For instance, we find a neutral 
fraction of ^ 37% for uh = 10"^ cm~^ and A^(H°) = 2.4 x 10^° cm"^. Of course, physically 
this would be imply a very large system with dimensions exceeding 100 kpc. 

6.2.2. CLOUDY simulations 

We have performed CLOUDY photoionization simulations with the aim of using 
the relative abundances of different ionization states of the same element to estimate 
the degree of ionization of our system. The analysis parallels the treatment presented 
by Lu et al. (1995). For the input radiation we used an attenuated power law ionizing 
spectrum computed by Madau (1992). The spectrum was normalized to a mean intensity of 
Jy = .195 X 10~^^ ergs s~^ cm~^ Hz~^ sr^^ at 1 Rydberg calculated at an average redshift of 
2;=2.46. The calculations assume a plane-parallel geometry with only one face illuminated 
and with incident radiation perpendicular to the surface, in contrast to the more realistic 
assumptions used to calculate H°/H in the previous section. The metallicity is fixed at 
[Z/H] = —0.5, the hydrogen volume density is varied, and the program terminates at a 
column density log[iV(H°)] = 20.38. Thus the neutral hydrogen column density is fixed at 
the observed value log[iV(H°)] = 20.38 while the total hydrogen column density is allowed 
to vary. Although this simulation is overly simplified (e.g. single side illumination, normal 
incidence, no dust depletion model), it does provide a good estimate of the degree of 
ionization. Figure |T5| presents the results of the simulations. Note that no conclusions 
concerning the degree of ionization can be drawn by comparing the relative abundances of 
the low-ions because their column densities track A^(H°) which is held constant. 

Given the CLOUDY results, we can use the observed Si^ to Si^'^ and Al"*" to Al^^ 
ratios to determine the ionization level of the gas. Since the AlII profile is heavily saturated, 
we take N{v) = 12.00 per pixel corresponding to Ia{v)/Ii{y) = 0.05 in the optical depth 
equation (Equation 1). Given the degree of saturation, this value yields a very conservative 
lower limit for [A1+/A1++] = log[A^(Al+)] - log[A^(Al++)]. Figure [T^ is a plot of (a) 
[Si"'"/Si'^^] and (b) [A1"'"/A1"'""'"] for 5 pixel bins over the entire low-ion region. The large dots 
represent the average value over the velocity regions defined by the vertical dashed lines and 
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the borders of the plot. Table 11 gives the values of [Si"'"/Si^"'"] and [A1"'"/A1'*"+] as well as 



the corresponding H^/H"*" ratio derived from the CLOUDY results for the velocity regions. 

With the exception of the region —40 km s~^ < i; < 40 km s^^ , we find [A1+/A1++] 
> 0.4 dex. This yields a neutral hydrogen fraction H^/H > 0.5. The velocity region v = 
[-100, 40] km s~^ corresponds to the strongest feature in all of the low-ion profiles and we 
expect the abundance of Al II to greatly exceed the value assumed above. The [Si+/Si^+] 
results are more difficult to interpret because the Si IV line-profile does not closely trace the 
low-ion or Al III profiles. According to the Sill to SilV ratio, the most highly ionized region 
is —218 km s^^ < v < —140 km s~^ which is in clear contradiction with the [Al^/Al^+j 
results. We believe that a significant portion of the SilV absorption is due to gas which 
is physically separate from the low-ion and A1++ gas. As such, we consider the [Si+/Si'^+] 
values to be lower limits. To summarize, our results suggest the hydrogen fraction of the 
gas 11*^/11 > 0.5 is conservative lower limit. Therefore, it is highly likely that the gas in the 
z = 2.46 damped system is mainly neutral. 

It is also likely that the neutrality of the gas is unaffected by the presence of internal 
sources of ionizing radiation such as O B stars. We checked this possibility by performing 
CLOUDY calculations with input from a black body with T = 4x10^ K. The resulting 
A1+/AU+ and Si^/Si"^^ ratios are indistinguishable from the ratios predicted by the 
external radiation field considered above. The results are similar because the sharp drop in 
ionizing fiux at the Lyman limit in the attenuated external spectrum (|Madau 1992|) mimics 



the exponential fall off of a black body. As a result the metal-line ratios indicate H is mainly 
neutral in the z = 2.46 absorber, for plausible sources of external and internal radiation. 



7. RESULTS 

This section presents the results from our abundance and kinematic analyses of the 
damped Lya system at z=2.462. We discuss the evidence for dust in this system and 
remark on the nearly constant abundances relative to Zn for three of the velocity features 
in our system. Finally, we compare and contrast the kinematics of the damped Lya system 
with several other systems, including another published HIRES damped Lya system (|Wolfe 
8t al. 19941). 
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7.1. Abundances of the ^=2.462 System 

Table |1^ lists the column density log[A^(X)], and the logarithmic abundance 



of element X relative to hydrogen normalized to solar abundances, [X/H] 
= log[A^(X)/A^(H)] - log[A^(X)/iV(H)]0, for the damped Lya system at 2=2.462. 
The abundance is derived assuming log[A^(H)] = 20.38 ± 0.045 (10% error) and standard 
solar abundances ( [Anders fc Grevesse 19891) . Our results indicate a relatively metal-rich 
system. If the three Zn+ features (discussed in §5.1) comprise 50% of the total Zn"*" 
abundance (an analysis of these components in the other low-ions suggests ~ 45%), we 
find [Zn/H] = —0.262. This is the most metal-rich of any damped Lya system for which 
accurate abundances have been determined at redshift z > 2.0. Furthermore, it has a higher 
metallicity (assuming Zn predicts metallicity; see below) than all but one of the systems 
measured in the most extensive survey of metallicity in damped Lya systems carried out so 
far (IFettmi et al. 1994|). 



Figure |T8| and Table |T3| present the abundances of the low-ions relative to Zn (assuming 
cosmic abundances) in the three velocity features where we could reliably measure the Zn 
abundance. The error bars are not entirely accurate (they are derived with the apparent 
optical depth method), but do serve as valuable guides. The overall variation of the 
abundances relative to Zn is very small for the 3 features and all of the elements are 
depleted with respect to Zn. The values of [X/Zn], however, do show a slight increase 
from the first to the third feature. We believe the most likely explanation lies in our 
measurements of the Zn abundance, particularly since features 1 and 2 are more likely to 
be blended with absorption from Cr II 2062 and because the trend is relatively systematic. 
Even given these minor differences, the relative abundances between Fe, Ni, Si, Cr and Zn 
are essentially constant over the 3 features (corresponding to nearly 200 km s~^ in velocity 
space). Therefore, we see little evidence for gas-phase abundance variations throughout our 
system. This observation indicates damped Lya systems are chemically well-mixed which 
further suggests they are detected at ages large compared to their internal dynamical time 
scales (i.e. the rotation period). 



7.2. Dust Depletion 

In this section we analyze variations of the gas-phase element abundances in damped 
Lya systems with (a) velocity and (b) condensation temperature. Such variations have 
been observed in the ISM of the Galaxy and have been used to infer properties of the dust 
responsible for element depletion. We wish to see whether similar effects are present at high 
redshifts. 
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7.2.1. Abundance Variations 

The absence of variations in [X/Zn] with respect to velocity in the z = 2.462 system 
is contrary to the presence of such variations in the ISM. Spitzer & Fitzpatrick (1993) and 
Spitzer & Fitzpatrick (1995) recently used GHRS spectra to infer these variations for the 
sightlines to the Galactic stars HD 93521 and HD 149881. These sightlines are relevant 
because the measured A^(H°) are similar to the value inferred for the damped system. 
Moreover, the depletion levels are similar, provided one interprets negative values of [X/Zn] 
in damped Lya systems to result from grain depletion. From Figure [l^ and Table [T^ we 
infer an upper limit of ~ 0.3 dex for the variation of [Si/Zn], [Fe/Zn], [Ni/Zn], or [Cr/Zn] 
across the 3 velocity features in the z = 2.46 damped system, while the variation for [Fe/Zn] 
in HD 93521 is ^ 0.9 dex and for [Cr/Zn] in HD 149881 is 0.6 dex ( [Spitzer and l^'itzpatrick| 



1993| , [Spitzer and I'ltzpatrick 1995] ). Fitzpatrick and Spitzer attribute the changes in [X/Zn] 



to the increase in dust destruction for clouds with increasing random velocity with respect 
to galactic rotation speed. 

Variations in [X/Zn] could also be due to density variations along the line of sight. 
Studies of the ISM ( [Jenkins 1987|) have shown correlations between the average Hydrogen 
volume density uh and the degree of depletion D relative to cosmic abundances for a 
variety of elements, including Si, Cr, and Fe. Jenkins fitted the ISM data with the following 
depletion curve: 



D = do + m[\ognH + 0.5] (2) 

where D is the logarithmic depletion relative to cosmic abundances (i.e., [X/H]), m is 
the slope, do is the value of D at log uh = —0.5 and hh is the Hydrogen volume density 
averaged over the line of sight to a given star. It is straightforward to show that the 
difference of two depletion values corresponding to two measurements {AD = D2 — Di) can 
be related to the ratio of the uh values of each measurement: 



As stated above the variation of Si, Cr, and Fe relative to Zn over the 3 features is 
only ^ 0.3 dex, where the largest variation is between the first and third features (the 
reddest and bluest). As suggested in §7.1 this variation is almost certainly systematic, 
most likely a result of blending with Cr II 2062. Therefore, the value can be considered an 
upper limit to the actual abundance variation over these 3 components. Table |I^ lists the 
depletion variations, the ISM slopes cited in Jenkins (1987), the corresponding predicted uh 
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variation, and la errors for Fe, Si, and Cr. We note the Cr measurements place the tightest 
constraints on the predicted variation of tih- Of course, the errors are rather large, but we 
can confidently predict that the variations of uh are no larger than ^ 1 dex provided this 
technique is apphcable. 

What suppresses variations of [X/Zn] in the damped system? If dust is present, then 
the sightline through the damped system must encounter an ISM in which nn varies by 
less than a factor of 10 along the line of sight. Moreover, the efficiency of grain destruction 
by shocks must be less efficient than in the ISM. Because the frequency of supernovae 
explosions account for the inhomogeneous nature of the ISM ( |McKee fc Ostriker 1977] ) and 
because supernovae are the main contributor of shocks in the ISM, then a lower frequency 
of supernova explosions could account for the lack of variations in the damped system. 
However, in the ISM a sightline of a few kpc typically encounters 4 orders of magnitude 
variation in uh (|Kulkarni fc Heiles 1987 ). A drastic reduction in the rate of supernova 
explosions would be required to reduce this to one order of magnitude. But models for the 
evolution of disk galaxies suggest that supernova input in the past was stronger not weaker 
than in the present (|Bruzual 1992|) . Other, nucleosynthetic considerations suggest that 
supernovae rates in the past were as least as frequent as they are at present ([Truran 1995|) . 
Therefore, while we cannot rule out this explanation altogether, it seems implausible. 

Therefore, the null detection of variations in [X/Zn] requires a weaker coupling of the 
depletion rate to the gas density, and a lower grain destruction rate by interstellar shocks 
than in the ISM. As a result, dust in the damped system is either absent or has significantly 
different properties from dust in the ISM. 



7.2.2. Condensation Temperature 

Figure |19| plots the gas phase abundances versus condensation temperature, Tc, for 
the observed low-ions from the damped Lya system (solid dots) and from the line-of-sight 
toward the Galactic star ( Oph ( |Uardelli 1994| ; open square^). The idea is to determine 



whether the anti-correlation between [X/H] and Tc observed in the ISM ( [Jenkins 1987| ) is 



also present in high-z damped Lya systems. Because Tq is the temperature at which half 
the gas phase atoms in a stellar atmosphere condense to solid form, the ISM anti-correlation 
is taken as evidence for dust formation in stellar atmospheres ( [b'ield 1974] ). Comparison 
between the ISM and damped diagrams shows some similarities. Specifically, the relative 
rankings of [Cr/H], [Si/H], [Fe/H], and [Ni/H] are about the same. However, the absolute 



■^The error bars for the ISM data are on the order of the size of the squares or smaller. 
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values are significantly higher in the damped system. In addition [Zn/H] in the ISM greatly 
exceeds all the other [X/H], but in the damped system [Zn/H] does not greatly exceed 
[Cr/H], [Fe/H], and [Ni/H], and may in fact be comparable to [Si/H]. The differences 
between the ISM and damped [X/H] may indicate the presence of gas in which the 
metallicity (inferred from undepleted [Zn/H]) is > —0.5 and has a lower dust-to-gas ratio 
to account for the smaller difference between [Zn/H] and [X/H] (cf. Pettini et al. 199"^ . 

On the other hand, Lu et al. (1995, 1996a, 1996b) interpret these patterns in terms 
of nucleosynthetic yields from type II supernovae. They point out that the damped Lya 
abundance patterns of N/0, Si/Fe, Cr/Fe, and Mn/Fe are consistent with the abundance 
patterns observed for population II halo stars which have been primarily enriched by 
type II supernovae ( [Wheeler et al. 1989| ). They also stress that Mn/Fe for damped Lya 
systems and halo stars disagrees with Mn/Fe inferred for the ISM, which is influenced by 
grain depletion. The only difficulties with the Lu et al. hypothesis stem from Zn/Fe. The 
quantity [Zn/Fe] > in damped Lya systems, while [Zn/Fe] ~ for stars with —3 < [Fe/H] 
< ( [Sneden et aLlMll ). 

Theorists working on chemical evolution have long been puzzled as to why [Zn/Fe] is 
independent of [Fe/H]. The closed box model by Malaney and Chaboyer (1996) predicts 
[Zn/Fe] to increase with [Fe/H] owing to the metallicity- dependent yield computed for 
type II supernovae explosions. These authors suggest that the constant [Zn/Fe] could 
be peculiar to the chemical evolution of the Galaxy, due to non-LTE effects in stellar 
atmospheres, or errors in the calculated yields. In addition, while the Malaney-Chaboyer 
model predicts [Zn/Fe] < 0, Hoffman et al. (1996) suggest neutrino driven winds in type II 
supernovae might result in [Zn/Fe] > 0. This effect may help explain the overabundance of 
Zn relative to Fe in damped Lya systems. Therefore, while some questions remain to be 
addressed, nucleosynthetic yields and grain depletion are equally plausible explanations for 
the abundance patterns observed in damped Lya systems. 

These competing explanations lead to at least two different ways for interpreting the 
metallicity of damped Lya systems. First, if the Lu et al. interpretation is correct and Zn 
is somehow overproduced with respect to Fe, than the metallicity of damped Lya systems 
is significantly lower (^0.7 dex) than estimated from the Zn abundances. In terms of the 
system analyzed in this paper, the metallicity would be —0.83 dex which would still be the 
highest metallicity observed in a damped Lya system with z > 2.0. On the other hand, if 
we accept depletion, the metalhcity of these systems is at least as high as the Zn abundance, 
and the Malaney-Chaboyer results indicate it could be even higher (possibly requiring an 
even greater level of depletion) . 
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7.3. Kinematics of the Damped Lya System (z=2A62) 

In this subsection we discuss the kinematics of the damped Lya system at 2;=2.462. We 
stress again the observed differences between the high and low-ion profiles and intercompare 
the velocity profiles of the 2;=2.462 system and several other damped systems. 

7.3.1. Comparison of the Low and High Ion Profiles 

Although all of the transitions from the z = 2.462 system span approximately the same 
velocity interval (~ 200 km s~^ ), there are marked differences between the low and high-ion 
line profiles. Similar to the damped Lya system toward PHL 957 discussed by Wolfe et al. 
(1994), the velocity profiles of all of the low-ions are visibly asymmetric with the component 
at the red edge being the strongest of the profile. In this case, the strongest component is at 
the red edge, where in PHL957 the strongest component is at the blue edge. Figure [2^ and 
Table |l^ give the column density by percent of total in 5 binned intervals in velocity space 
corresponding to visible features in the line profiles. This serves as quantitative evidence of 
the asymmetry in the low-ions and highlights the inherent differences between the low and 
high-ion profiles. The overall profile of the high-ions has nearly the opposite shape of the 
low-ion profiles with the strongest feature on the blue edge. In addition, the profile of the 
high-ions is smoother with less structure than the low ions. 

We contend these differences, particularly those evident in the Si II and Si IV profiles, 
indicate the absorption associated with the high-ions does not originate in the same region 
as that associated with the low-ions. The most compelling evidence for this hypothesis lies 
in the extreme difficulty encountered in obtaining an accurate fit to the Si IV profile by 
altering the column densities of the calculated Si II components while holding the redshift 
and b values constant. Both the smoothness and the overall shape of the Si IV profile 
require significantly different b values and argue against the inclusion of so many thin, 
resolved components. 

7.3.2. The z =2.457 subsystem 

The neighboring system at z=2.457 was unidentified in the lower resolution data. 
It has similar line profile structure as the 2;=2.462 system with a large velocity interval, 
smoother high- ion profiles and low- ion profiles with more velocity structure. It does not, 
however, exhibit the edge leading asymmetry apparent in the z=2.462 system. In fact, the 
low- ions appear relatively symmetric about f = 70 km s~^ . We contend these differences in 
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the kinematic characteristics result from the fact that the Z—2A57 system is not a damped 
Lya system but most hkely an ionized Lyman hmit system in which the gas kinematics are 
not determined by rotation. 

7.3.3. Comparison with PHL 957 

Contrasting the damped Lya metal transitions with those toward PHL 957 aX z — 
2.309, we note several important differences. First, the velocity interval of the Q0201+365 
damped system is significantly larger (200 km s~^ vs. 50 km s~^ ). Secondly, the low- ions 
of the Q0201-I-365 system exhibit more velocity structure, an effect which is not due to 
differences in resolution. Finally, the asymmetric shape is not as prevalent in our damped 
system possibly because of its larger velocity interval. That is, the gas is more evenly 
distributed in Q020 1+365 than in the PHL957 absorber. In terms of the thick rotating 
disk model, these differences may be explained by differences in the inclination angle of the 
line-of-sight with respect to the disk. 

7.3.4- Comparison with System at z — 2.325 

A comparison of the damped Lya kinematics with the absorption system aX z = 2.325 
toward Q0201+365 further stresses the damped Lya characteristics. The z = 2.325 system 
(a possible Lyman limit system) also spans ~ 200 km s"-*^ in velocity space, but does 
not possess the same edge-leading asymmetric profile observed in the damped systems. 
Furthermore, the high- ion profiles at 2;=2.325 trace the low-ion profiles more closely than 
in the damped Lya systems. For instance, Al III was successfully tied to the low-ion 
VPFIT solution. These characteristics are all consistent with an explanation depicting the 
2;=2.325 system as a multi-cloud system in which the velocities are random rather than 
systematic as in the case of rotation. Moreover the gas is more evenly distributed. In short, 
the z = 2.325 system more closely resembles the ^=2.457 system which is also a likely 
Lyman limit system. Perhaps we are observing kinematic characteristics which differentiate 
damped Lya from Lyman limit systems. However, the statistics of small numbers makes 
this observation a speculation rather than a conclusion. 



-22- 



8. CONCLUSIONS 

This paper presented HIRES spectra obtained with the Keck 10m telescope of 
absorbing gas toward toward Q0201+365. We identified over 80% of the absorption features 
and have analyzed several of the more interesting metal-line systems. We have focused 
on the damped Lya system at z=2.462 as part of an ongoing program to investigate the 
chemical content and kinematics of damped systems within the redshift interval z ^ 2 — A. 
We summarize our results as follows. 

(1) Based on the analysis of ionization simulations, we predict the damped Lya system 
to be significantly neutral. Although it is possible the system is partially ionized, our 
analysis predicts the metals are all essentially in the singly ionized state, and that the total 
hydrogen column density is well within a factor of two of the adopted value from the A^(H°) 
measurement. 

(2) A hidden component analysis of the Ni II 1741, 1751 transitions did not reveal any 
significant hidden saturated components. We expect this to hold true for the other low-ion 
transitions. 

(3) With the VPFIT least squares line profile fitting package we have measured ionic 
column densities for the damped Lya system at 2;=2.462 (as well as several other systems). 
We performed similar measurements with the apparent optical depth method and found the 
two results to be in agreement, further eliminating the possibility of hidden line saturation. 

(4) We measured the following abundances of Si, Fe, Cr, and Ni for the damped Lya 
system: [Si/H] = -0.376 ± 0.052, [Fe/H] = -0.830 ± 0.051, [Cr/H] = -0.902 ± 0.064, and 
[Ni/H] = —1.002 ± 0.054. We placed limiting values on the abundances of the s-process 
elements Pb (3ct detection) and Ge (upper limit), [Pb/H] = 2.233 ± 0.121 and [Ge/H] < 
0.664, and a lower limit value on the abundance of Zn, [Zn/H] > —0.562. Based on the 
VPFIT solution of the low-ions, we expect the metallicity is [Z/H] ^ —0.262. This damped 
Lya system has the highest metallicity measured to date at 2; > 2.0. 

(5) Comparing individual features of the damped Lya system, we find the relative 
abundance between Si, Fe, Cr, Ni and Zn remains nearly constant throughout our system 
(Figure |l^). This suggests a well mixed system with an age large compared to the internal 
dynamical time scale at the epoch of detection. 

(6) We have used the relatively minor variations observed in the Si, Cr, and Fe 
abundances relative to Zn to place limits on the expected variation in the Hydrogen volume 
density throughout the damped Lya system, having assumed the presence of dust grains 
and the Jenkins relation ( [Jenkins 19871 ). Our measurements of [Cr/Zn] place a maximum 
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variation of hh at ^ 1 dex. The lack of uh (and [X/Zn]) variations could be evidence of 
weaker supernova input in the past, but we believe they are more likely due to the absence 
of grains with the properties of dust found in the ISM. 



(7) Plotting the measured abundances versus condensation temperature (Figure |T9|), 
we do find evidence for a depletion pattern, but the overall depletion level of Si, Fe, Cr and 
Ni with respect to Zn is indicative of a relatively dust free ISM cloud. Although gas with a 
lower dust-to-gas ratio than evident in the ISM can account for the pattern, one can also 



explain the pattern in terms of nucleosynthetic yields from type II supernovae (|Lu, Sargent 



|fc Barlow 1995| ). Both explanations are problematic and are under continued debate. 



Determining the proper explanation is particularly important as they predict different 
metallicities for damped Lya systems, which will significantly affect the investigation of 
galactic chemical evolution. 

(8) The low-ion profiles of the damped system exhibit an edge leading asymmetry 
as predicted by a simple model of rotation. The shape is similar to the other damped 
system observed with HIRES (PHL 957; [Wolfe et al. 19941) , though the velocity interval is 
significantly greater. 
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UT Date Exposure Wavelength Resolution 
Time (s) Coverage ( A) (km ) 



1994 Sep 15 9600 
1994 Sep 30 11950 
1994 Oct 1 13030 



4720 - 7130 7.2 - 8.0 
4790 - 7180 7.0 - 8.0 
4790 - 7180 7.2 - 8.1 



Table 1: JOURNAL OF OBSERVATIONS 
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4876.08 


0. 


.0190 





.0028 










4895.98 


0, 


(1982 





.0087 


CI(1656) 


1 


.955 


72 




















4924.10 


0, 


.1703 


0. 


.0073 


CII(1334) 


2. 


.690 




4938 




ga 






A1II(1671) 


1 


.955 




4940 




ga 






Mgl(2852) 


0. 


.732 




4951.84 


0, 


.2428 


0. 


.0112 


SiIV(1393) 


2. 


.553 




4984.16 


0, 


.2462 


0. 


.0152 


SiIV(1402) 


2. 


.553 


71 




















5005.84 


1, 


.3215 


0. 


.0112 


CIV(1548) 


2. 


.233 




5014.00 


0, 


.8942 


0. 


.0101 


CIV(1550) 


2. 


.233 


70 




















5076.49 


1, 


.3223 


0. 


.0153 


SiII(1526) 


2. 


.325 


69 




















5139.61 


0, 


.3223 





.0059 


CIV(1548) 


2 


.320 




5142.70 


0, 


.7369 


0. 


.0094 


SiIV(1393) 


2. 


.690 




5148 




ga 






CIV{1548) 


2. 


.325 




5148.15 




ga 






CIV(1550) 


2. 


.319 




5156.16 


2. 


.0353 


0. 


.0133 


CIV(1550) 


2. 


.325 




5172.86 




Sb 






Sky Emiss 








5175.62 


0, 


.8012 





.0136 


SiIV(1402) 


2 


.690 




5177.77 


0, 


.0241 





.0035 










5189.73 


0, 


.0407 


0. 


.0054 


CI(1560) 


2. 


.325 




5200.18 


0, 


.0301 


0. 


.0054 


Sky?? 








5201.43 


0, 


.0609 


0. 


.0057 


CIV(1548) 


2. 


.360 
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Table 2 — Continued 



Order 


A 

(A) 




w 
(A) 




ID 




ofcs 




'=191 n 07 


0. 


.0330 




0047 




2 


.360 




KnOA Oft 




.0337 


Q 


0049 


oKy ADS 








( U.o4 







Q 


0076 










5272.51 




2489 
















Q 


1389 





.0077 














6663 





.0177 




2 


457 


67 




















5285.95 


3. 


.3202 





.0121 


OITT/'i co^;^ 
01iH10.iD } 




.462 






2 


.4749 






V ^lo4c4 ) 


2 


.423 








7411 




OlM 


V (loou } 




.423 




OO^D.UO 


0. 


.3089 


0, 


.0069 


oilV (^loyo ) 


2 


.821 






° 




0, 


.0131 


q;TT/i sJOSJ"! 








oo4 . 


0. 


2201 


0, 


.0148 


T.'.iTT ^ 1 AO\ 

i^eii(^iuUo } 




one 
.325 




5352.75 


3, 


.2168 


0, 


.0166 


CIV(1548) 


2, 


.457 




5357.45 


2, 


.3480 


0, 


.0169 


CIV(154S) 


2, 


.462 


66 




















5367.53 


2 


.1787 




0194 


CIV{1550) 


2 


.462 




DO / O.UO 






0, 


.0091 


K^L V ±oou ) : 








OOOO.U4 


n 




0, 


.0130 


^1. V ±040 ) 




477 




5391.61 




1891 





.0124 


^1. V ^^ IDOU ) 




477 




6419.46 


0, 


.4062 





.0097 


CIV (1548) 


2, 


.500 




6426.90 










Sky Abs 








5428.50 





.2452 





,0112 


CIV(1550) 


2 


.500 


65 




















o4oo. la 


0. 


.0286 





.0039 










6481.60 


1, 


.3494 





.0128 


A1III(1854) 


1 


.955 




6600.98 


1, 


.6163 





.0100 


CIV{1548) 


2 


.553 




6606.02 


0. 


.7622 





.0135 


A1III(1862) 


1 


.955 




6610.20 


0. 


.9414 





.0119 


CIV(1550) 


2 


.553 




6612.77 


0. 


.0255 





.0042 








64 




















6666.32 


1. 


.3290 





.0144 


Alll^lD 11.) 


2 


.325 




6660.62 


0. 


.2871 





.0130 


FeII(1608) 


2 


.457 




6668.31 


2. 


.6714 





.0117 


rell(^iDOo ) 


2 


.462 




oo7y.22 




ga 






remiDll ) 


2 


.462 




5607.60 


0. 


.0793 


0, 


.0079 








63 


6622.07 


0, 


.0397 


0, 


.0062 








62 


6713.13 


1, 


.0943 


0, 


.0136 


CIV{1548) 


2, 


.690 




6722.31 


0, 


.6010 


0, 


.0125 


CIV(1550) 


2, 


.690 




6776.46 


2, 


.1480 


0, 


.0137 


A1II(1671) 


2, 


.457 


61 


6783.40 


6, 


.0321 


0, 


.0164 


A1II(1671) 


2, 


.462 




5804.61 





,5622 





,0091 


FeII(2344) 


1 


.476 




6846.16 


0, 


,0226 


0, 


.0044 










6872.81 


0. 


.0190 


0, 


.0037 


FeII(2374) 


1 


.473 




5874.03 


0, 


.1010 


0, 


.0072 


FeII(2374) 


1 


.474 




5879.60 


0, 


,2519 





.0072 


FoII(2374) 


1 


.476 


60 




















5901.21 


0, 


,4404 


0, 


.0099 


FoII(2382) 


1 


.476 




5916.32 




ga 






CIV{1548)? 


2 


.821 




5919 




ga 






NiII(1709) 


2, 


.462 




6076.31 


0, 


.0191 


0, 


.0033 


CrII(2066) 


1 


.955 


69 




















6010.37 


0, 


,0947 


0, 


.0127 


SiII(1808) 


2, 


.325 




6028.97 


0, 


,3928 


0, 


.0113 


NiII(1741) 


2, 


.462 




6043.72 










Ink Blot 








6064.91 


0, 


.2400 


0, 


.0110 


NiII(1761) 


2, 


.462 
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Table 2 — Continued 



Order A 


Vv 




in 


^abB 




(A) 










58 














6166.72 


0.6211 


0.0140 


A1III(1854) 


2.325 


57 














6259.01 


0.6740 


0.0140 


SiII(1808) 


2.462 




6280 - 






Sky Abs 




56 
















0.3591 


0.0090 


FeII(2344) 


1.699 




6331.86 






Sky Abs 






6366 






Sky Abs 






6380.61 






Sky Abs 






ooyo.oo 


0.2889 


0.0095 


FcII{2586) 


1.473 




6400 '^^ 


0.0246 


0.0040 


FeII{2586) 


1.474 




6403.12 


0.0902 


0.0061 


FeII(2586) 


1.475 




6404.76 


0.5519 


0.0089 


FeII(2586) 


1.476 




6411.14 


0.2570 


0.0113 


FbII(2374) 


1.699 


55 














R4^8 1 


0.9346 


0.0102 


FcII(2600) 


1.476 




6440.63 


0.2097 


0.0109 


AlIII(18(i3) 


2.457 




6448.40 


0.9757 


0.0148 


A1III{1863) 


2.462 




6470 - 






Sky Abs 




54 












53 












52 














6813.84 


0381 


0.0053 








6869 - 






Band Abs 




51 














6942.03 


0.9711 


0.0084 


MgII(2803) 


1.476 




6944.50 






Sky Abs 






6946.10 






Sky Abs 






6958.54 






Sky Abs 






6963.62 






Sky Abs 






6982.78 


0.2482 


0.0105 


FbII(2586) 


1.699 




6987.39 






Sky Abs 






6988.89 






Sky Abs 






6991.29 






Sky Abs 






6995.81 






Sky Abs 






7001.26 


s" 




Sky Abs 






7007.03 






Sky Abs 






7016 






FeII{2374) 


1.955 




7018.78 






FeII(2600) 


1.699 




7039.67 


1.6433 


0.0080 


FeII(23S2) 


1.955 


50 














7118.18 


0.2471 


0.0159 


CrII{2056) 


2.462 




7140.5 






CrII{2062) 


2.462 




7140.5 


B» 




ZnII(2062) 


2.462 




7144 - 






H2O Abs 




a 


Line blend 


ing 








b 


Night-sky absorption or emission 



features 
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Table 3. METAL TRANSITIONS 



iransixion 


\ ( °A\ 
^rest {^) 


J 


Nill 1370 


1370.132 


0.1309 


SilV 1393 


1393.755 


0.528 


SnII 1400 


1400.400 


0.71 


SilV 1402 


1402.770 


0.262 


Gall 1414 


1414.402 


1.8 


Sill 1526 


1526.707 


0.2303 


CIV 1548 


1548.195 


0.1908 


CIV 1550 


1550.770 


0.09522 


Cell 1602 


1602.4863 


0.135 


Fell 1608 


1608.4449 


0.05545 


AlII 1670 


1670.7874 


1.88 


Pbll 1682 


1682.15 


0.156 


Nill 1741 


1741.549 


0.1035 


Nill 1751 


1751.910 


0.0638 


Sill 1808 


1808.0126 


0.00218 


Aim 1854 


1854.716 


0.539 


Aim 1862 


1862.790 


0.268 


Znll 2026 


2026.136 


0.515 


CrII 2056 


2056.254 


0.1403 


CrII 2062 


2062.234 


0.1049 


Znll 2062 


2062.664 


0.2529 


Fell 2344 


2344.214 


0.1097 


Fell 2374 


2374.4612 


0.02818 


Fell 2382 


2382.765 


0.3006 


Fell 2600 


2600.1729 


0.2239 


Mgll 2796 


2796.352 


0.6123 


Mgll 2803 


2803.531 


0.3054 
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Table 4a. FIT FOR ^=2.462 - LOW IONS 



Comp 


z 


X 10" = 


b 

(km ) 


(km ) 


ION 


log N 
(cm-2) 


(cm- = ) 


1 


2.460470 


2.0 


4.59 


2.19 


Fe+ 


13.11 


0.15 


2 


2.460598 


0.3 


4.44 


0.47 


Fe+ 


14.11 


0.05 












Si+ 


14.41 


0.04 












Ni+ 


12.56 


0.04 












Cr+ 


12.02 


0.10 


3 


2.460807 


0.5 


5.91 


1.01 


Fe+ 


13.99 


0.13 












Si+ 


14.36 


0.16 












Ni+ 


12.30 


0.32 












Cr + 


12.18 


0.12 


4 


2.460944 


5.9 


12.47 


4.97 


Fe + 


13.80 


0.24 












Si+ 


14.21 


0.25 












Ni+ 


12.54 


0.22 












Cr+ 


11.61 


0.55 


5 


2.461347 


1.4 


13.94 


1.17 


Fe+ 


13.99 


0.04 












Si+ 


14.52 


0.06 












Ni+ 


12.64 


0.07 












Cr + 


12.28 


0.12 


6 


2.461426 


0.5 


3.00 


1.22 


Fe + 


13.74 


0.10 












Si+ 


13.94 


0.17 












Ni + 


12.00 


0.19 












Cr + 


11.69 


0.31 




^.^Ol 1 oO 








Fe+ 
















Si+ 


14.76 


0.03 












Ni+ 


12.85 


0.03 












Cr+ 


12.23 


0.08 


8 


2.461905 


0.4 


3.24 


0.89 


Fe+ 


13.92 


0.08 












Si+ 


14.33 


0.07 












Ni+ 


12.45 


0.07 












Cr+ 


11.96 


0.13 


9 


2.462091 


1.4 


8.23 


2.59 


Fe+ 


13.68 


0.11 












Si+ 


14.19 


0.14 












Ni+ 


12.30 


0.12 












Cr+ 


11.06 


1.15 


10 


2.462266 


1.5 


6.22 


2.42 


Fe+ 


13.51 


0.19 












Si+ 


14.13 


0.19 












Ni+ 


11.94 


0.30 












Cr+ 


11.02 


1.22 


11 


2.462494 


2.2 


11.31 


2.43 


Fe + 


13.99 


0.10 












Si+ 


14.59 


0.11 












Ni+ 


12.59 


0.11 












Cr + 


12.27 


0.10 


12 


2.462594 


0.6 


3.00 


1.49 


Fe+ 


13.88 


0.18 












Si+ 


14.17 


0.17 












Ni+ 


12.18 


0.18 


13 


2.462818 


0.6 


14.58 


0.56 


Fe+ 


14.43 


0.02 












Si+ 


14.92 


0.02 












Ni+ 


13.00 


0.02 












Cr+ 


12.57 


0.05 
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Comp z (7z 

( X 10-5) 



u 




Ion 


lU^ i V 




(km s^^ ) 


(km s^^ ) 




( cm"^) 


( cm"^ 


18.33 


3.03 




11.72 


0.53 








13.47 


0.07 


43.14 


8.77 




12.81 


0.12 








13.95 


0.04 


9.43 


0.99 




12.77 


0.07 








13.45 


0.10 


16.65 


1.50 




13.45 


0.04 








14.19 


0.13 


11.63 


0.64 




13.54 


0.04 








13.10 


0.24 


82.54 


22.64 




13.36 


0.14 






C+3 


12.85 


0.38 


15.01 


1.61 




12.81 


0.08 








13.41 


0.48 


11.01 


5.97 




11.67 


0.52 


38.26 


11.23 




12.81 


0.31 






C+3 


13.35 


0.11 


5.38 


5.93 




11.59 


0.34 






C+3 


13.63 


0.12 



1 2.459346 2.5 

2 2.459713 4.4 

3 2.460433 1.0 

4 2.460753 1.2 

5 2.461035 0.8 

6 2.461647 22.3 

7 2.461863 0.9 

8 2.462229 3.5 

9 2.462762 9.1 

10 2.462908 2.1 



Table 4b: FIT FOR ^=2.462 - HIGH IONS 
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Table 5a. FIT FOR z=2.457 - LOW IONS 







/T 

^ Z 

/\ J.VJ 


h 

U 

1 JVill O I 


1 JVlll o 1 


ION 


iUg, iV 

cm ) 


*^logiV 

( f*TY\ ^ 


1 


2.456041 


0.5 


3.95 


0.80 


Si+ 


12.53 


0.06 












Fe+ 


12.34 


0.20 












A1+ 


11.81 


0.04 


2 


2.456295 


0.3 


9.58 


0.49 


Si+ 


13.40 


0.02 












Fe+ 


13.19 


0.04 












A1+ 


12.60 


0.02 


3 


2.456540 


1.5 


5.44 


2.57 


Si+ 


12.54 


0.15 












Fe+ 


12.30 


0.29 












A1+ 


11.78 


0.12 


4 


2.456679 


0.6 


3.30 


1.30 


Si+ 


13.03 


0.05 












Fe+ 


12.78 


0.12 












A1+ 


12.19 


0.06 


5 


2.456844 


1.6 


8.95 


2.88 


Si+ 


12.98 


0.18 












Fe+ 


12.95 


0.14 












A1+ 


12.26 


0.16 


6 


2.457079 


1.5 


12.83 


1.78 


Si+ 


13.41 


0.06 












Fe+ 


13.04 


0.10 












A1+ 


12.65 


0.06 


7 


2.457336 


0.6 


3.00 


1.15 


Si+ 


13.21 


0.11 












Fe+ 


12.78 


0.09 












A1+ 


12.28 


0.11 


8 


2.457460 


0.6 


3.00 


1.22 


Si+ 


13.16 


0.10 












Fe+ 


12.63 


0.12 












A1+ 


12.22 


0.09 


9 


2.457689 


6.4 


25.16 


7.00 


Si+ 


13.04 


0.13 












A1+ 


12.37 


0.11 


10 


2.457858 


0.3 


6.37 


0.57 


Si+ 


13.35 


0.04 












Fe+ 


13.11 


0.05 












A1+ 


12.42 


0.05 


11 


2.458100 


1.1 


4.65 


2.29 


Si+ 


12.40 


0.15 












Fe+ 


12.17 


0.30 












A1+ 


11.42 


0.25 


12 


2.458406 


0.8 


14.96 


1.20 


Si+ 


12.80 


0.05 












A1+ 


12.09 


0.03 
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Comp 


z 


( X 10-->) 


b 

(km ) 


Oh 

(km ) 


Ion 


log 

(cm-2) 


(cm-^: 


1 


2.456343 


0.7 


24.36 


0.78 




12.98 


0.02 














13.86 


0.02 


2 


2.457273 


0.5 


34.58 


0.70 




13.65 


0.01 














14.54 


0.02 


3 


2.457873 


0.4 


8.77 


0.61 




12.98 


0.02 














13.75 


0.05 


4 


2.458106 


0.9 


7.53 


1.17 




12.45 


0.06 














13.32 


0.07 


5 


2.458454 


1.5 


30.67 


1.98 




12.67 


0.04 












C+3 


13.77 


0.03 



Table 5b: FIT FOR z=2.457 - HIGH IONS 
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Table 6a. FIT FOR z=2.325 - LOW IONS 



Comp 


z 


X 10~^ 


b 

(km ) 


(km ) 


ION 


log JV 
(cm-2) 


(cm-2) 


1 


2.323701 


0.6 


6.87 


0.99 


A1+ 


11.52 


0.06 














11.97 


0.05 












Fe+ 


11.87 


0.78 












Si+ 


12.73 


0.04 


2 


2.323896 


0.3 


3.39 


0.57 


A1+ 


13.00 


0.25 












Al+2 


12.59 


0.03 












Fe+ 


13.32 


0.04 












Si+ 


13.92 


0.24 


3 


2.324033 


0.3 


3.37 


0.69 


A1+ 


12.83 


0.23 












Al+2 


12.48 


0.03 












Fe+ 


13.30 


0.04 












Si+ 


13.81 


0.26 


4 


2.324173 


1.0 


3.41 


1 .86 


Al+= 


11.81 


0.07 


5 


2.324195 


2.8 


9.92 


3.16 


A1+ 


11.63 


0.11 












Si+ 


12.22 


0.17 


6 


2.324663 


0.2 


6.24 


0.26 


A1+ 


12.28 


0.02 












Al+2 


12.02 


0.04 












Fe+ 


12.71 


0.13 












Si+ 


13.21 


0.02 


7 


2.324884 


1.2 


2.80 


2.97 


A1+ 


11.11 


0.12 












Si+ 


11.98 


0.15 


8 


2.325078 


0.7 


3.12 


0.88 


A1+ 


11.97 


0.05 












Al+2 


11.92 


0.06 












Fe+ 


12.47 


0.22 












Si+ 


13.06 


0.06 


9 


2.325180 


0.4 


2.27 


0.78 


A1+ 


13.03 


0.70 












Al+2 


11.99 


0.06 












Fe+ 


13.36 


0.07 












Si+ 


14.30 


0.76 


10 


2.325560 


0.6 


5.66 


1.02 


A1+ 


11.45 


0.06 












Al+2 


11.65 


0.09 












Fe+ 


12.01 


0.56 












Si+ 


12.47 


0.06 


11 


2.325796 


1.5 


4.88 


2.55 


A1+ 


11.20 


0.10 


12 


2.326045 


0.3 


4.68 


0.39 


A1+ 


12.16 


0.03 












Al+2 


12.11 


0.04 












Fe+ 


12.46 


0.20 












Si+ 


13.34 


0.03 


13 


2.326198 


0.5 


6.88 


0.54 


A1+ 


12.08 


0.03 












Al+2 


12.09 


0.04 












Fe+ 


12.75 


0.12 












Si+ 


13.21 


0.03 
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Comp 


z 


( X 10-^) 


b 

(km ) 


(km ) 


Ion 


log iV 
( cm"^) 


( cm~^) 


1 


2.319744 


0.2 


9.11 


0.29 




13.64 


0.01 


2 


2.323508 


4.4 


15.48 


3.55 


C+3 


13.13 


0.13 


3 


2.323736 


1.2 


9.13 


1.83 


C+3 


13.47 


0.16 


4 


2.324048 


1.4 


24.68 


2.38 


C+3 


14.25 


0.04 


5 


2.324569 


2.4 


15.07 


3.74 


C+3 


13.56 


0.13 


6 


2.324825 


1.9 


11.15 


2.60 


C+3 


13.41 


0.18 


7 


2.325100 


3.8 


16.40 


4.30 


Q+3 


13.25 


0.12 


8 


2.325624 


0.7 


18.12 


1.23 




13.47 


0.02 


9 


2.326123 


0.2 


14.21 


0.25 


Q+3 


14.23 


0.01 



Table 6b: FIT FOR ^=2.325 - HIGH IONS 



)mp 


z 


( X 10-'^') 


b 

(km ) 


Oh 

(km s""*^ ) 


Ion 


log iV 
( cm~^) 


( cm~^ 


1 


1.475760 


1.5 


1.21 


5.88 


Fe+ 


12.24 


0.18 


2 


1.475858 


0.7 


4.27 


1.87 


Fe+ 


12.71 


0.07 


3 


1.476039 


0.4 


3.59 


0.52 


Fe+ 


13.65 


0.03 


4 


1.476115 


0.8 


2.44 


1.37 


Fe+ 


13.22 


0.06 


5 


1.476182 


0.9 


2.05 


1.86 


Fe+ 


12.92 


0.08 


6 


1.476301 


0.1 


2.70 


0.22 


Fe+ 


13.64 


0.04 



Table 7: FIT FOR z=1.476 
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Transition Apparent VPFIT Adopted 



Si IV 1393 


14, 


,039 ± 


0, 


,006 


14.050 ± 0.038 


14, 


,045 


±0, 


.038 


C IV 1550 


14. 


.617 ± 


0. 


,005 


14.615 ± 0.062 


14. 


,616 


±0, 


.062 


C I 1560 


13. 


,228 ± 


0. 


,059 




13. 


,228 


±0, 


.059 


Fe II 1608 


15. 


.033 ± 


0. 


,004 


15.086 ± 0.024 


15. 


,060 


±0, 


.024 


Al II 1671 




















Ni II 1709 


13. 


,785 ± 


0. 


,020 




13. 


,785 


±0, 


.020 


Ni II 1741 


13. 


,638 ± 


0. 


,013 


13.631 ± 0.030 


13. 


,628 


±0, 


.030 


Ni II 1751 


13. 


,610 ± 


0. 


,020 












Si II 1808 


15. 


,561 ± 


0. 


,010 


15.546 ±0.026 


15. 


,554 


±0, 


.026 


Al III 1862 


13. 


,606 ± 


0. 


,008 




13. 


,606 


±0, 


.008 


Cr II 2056 


13. 


,158 ± 


0. 


,030 


13.138 ±0.040 


13. 


,148 


±0, 


.040 


Cr II 2062 


13. 


,401 ± 


0. 


,025 




13. 


,401 


±0, 


.025 


Zn II 2062 


12. 


,468 ± 


0. 


,046 




12. 


,468 


±0, 


.046 



Pb II 1682 12.66 ±0.114 12.66 ±0.114 

Ge II 1602 < 12.65 < 12.65 



Table 8: IONIC COLUMN DENSITIES FOR z = 2.462 



Note. — Values reported in logarithmic space have deceptively small errors. For instance, the value for 
A''(Pb) is not a 5(7 detection and that for A^(Ge) is not even a la detection. 



Transition Apparent VPFIT Adopted 



Si IV 1393 
Si IV 1403 
C IV 1550 

Si II 1526 
Fe II 1608 
Al II 1671 
Al III 1862 



13.839 ±0.003 
13.851 ±0.007 
14.686 ± 0.041 

14.129 ±0.006 
13.839 ± 0.022 

13.321 ±0.004 
12.920 ±0.028 



13.85 ±0.01 

14.74 ±0.01 

14.18 ±0.02 
13.85 ± 0.03 
13.36 ±0.02 



13.85 ±0.01 
13.85 ±0.01 
14.737 ±0.010 

14.18 ±0.02 
13.843 ±0.018 

13.36 ±0.02 
12.920 ±0.028 



Table 9: IONIC COLUMN DENSITIES FOR z = 2.457 
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Loiiipoiiout 


logio A (HI) 


- abb 


7/1 

b (km s 


1 


18.38 


2.460470 


4.59 


2 


19.38 


2.460598 


4.44 


3 


19.26 


2.460807 


5.91 


4 


19.07 


2.460944 


12.47 


5 


19.26 


2.461347 


13.94 


6 


19.01 


2.461426 


3.00 


7 


19.47 


2.461735 


8.56 


8 


19.18 


2.461905 


3.24 


9 


18.95 


2.462091 


8.23 


10 


18.78 


2.462266 


6.22 


11 


19.26 


2.462494 


11.31 


12 


19.15 


2.462594 


3.00 


13 


19.70 


2.462818 


14.58 



14 


17.61 


2.456041 


3.95 


15 


18.46 


2.456295 


9.58 


16 


17.57 


2.456540 


5.44 


17 


18.05 


2.456679 


3.30 


18 


18.22 


2.456844 


8.95 


19 


18.31 


2.457079 


12.83 


20 


18.05 


2.457336 


3.00 


21 


17.90 


2.457460 


3.00 


22 


18.38 


2.457858 


6.37 


23 


17.44 


2.458100 


4.65 



Table 10: HI COMPONENTS IN Lya PROFILE (^=2.46) 
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Velocity (km s ^ ) 


[Si+/Si3+] 


[H7H+]si 


[A1+/A1++] 


[HVH+]ai 


[H7H+]ai 


-218 <v < -140 


1.0 


-0.16 


0.42 


-0.06 


-0.06 


-140 <v < -110 


1.7 


0.12 


0.47 


0.00 


0.12 


-110 <v < -40 


1.8 


0.13 


0.43 


-0.04 


0.13 


-40<v < 40 


2.1 


0.31 


0.27 


-0.30 


0.31 



Table 11: IONIZATION LIMITS 



Note. — All values are conservative lower limits 



Metal logioiV(X) (cm-2) [X/H] 

1^ 15.060 ±0.024 -0.830 ±0.051 

Ni 13.628 ±0.030 -1.002 ±0.054 

Al SATU 

Si 15.554 ±0.026 -0.376 ± 0.052 

Cr 13.158 ±0.030 -0.902 ± 0.064 

Zn > 12.468 ± 0.046 > -0.562 ± 0.064 

C 14.616 ±0.062 -2.324 ±0.077 

Si 14.045 ±0.038 -1.885 ±0.059 

Pb 12.66 ±0.114 2.23 ±0.121 

Ge < 12.65 < 0.644 



Table 12: ABUNDANCES FOR z = 2.462 



[X/Zn] Feature 1 Feature 2 Feature 3 

-144 ^ -110 km/s -40 ^ -9 km/s -9 ^ 35 km/s 

Fe -0.856 ±0.185 -0.679 ±0.175 -0.508 ±0.243 

Ni -0.939 ±0.204 -0.831 ±0.186 -0.638 ±0.249 

Si -0.388 ±0.197 -0.157 ±0.180 -0.035 ± 0.246 

Cr -0.833 ±0.261 -0.729 ±0.220 -0.534 ± 0.274 



Table 13: DEPLETION RELATIVE TO Zn FOR z = 2.462 



-38- 



Metal AD 

~~Fe 0.348 ±0.305 

Si 0.353 ±0.315 

Cr 0.299 ± 0.378 



0.38 ±0.05 8.2 ±9.7 
0.49 ±0.15 5.3 ±6.6 
0.50 ±0.11 4.0 ±5.3 



Table 14: uh VARIATIONS FOR z=2A62 



"m values are taken from Jenkins 1987 



Metal Comp 1 Comp 2 Comp 3 Comp 4 Comp 5 

-211^-181 -181^-144 -144 ^-110 -110^-40 -40 ^35 

Fe oil 014 oil 025 039 

Ni 0.09 0.14 0.13 0.27 0.38 

Si II 0.10 0.11 0.11 0.28 0.42 

Cr 0.11 0.15 0.14 0.23 0.40 

Si IV 0.17 0.46 0.09 0.14 0.10 



Table 15: PERCENT OF TOTAL ABUNDANCES FOR z = 2.462 
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Fig. 1.— Keck HIRES spectra of Q0201+365 at a resolution of 8 km s"^ and a SNR of 
^ 30. All 26 orders (only 2 are shown here) are identified by the echelle order intrinsic to 
the Keck HIRES spectrograph. The dotted line is a Icr error array derived assuming Poisson 
statistics. 
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Fig. 2. — Velocity profiles of those metal-line transitions from the system at z=2A62 with 
successful VPFIT solutions. The dashed vertical line corresponds to z=2A628. The marks 
above the Fe II 1608 and Si IV 1393 transitions indicate the velocity centroids for VPFIT 
solutions as listed in Table 4. In all velocity plot figures with VPFIT solutions, the leftmost 
mark is component 1. 
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Fig. 3. — Velocity profiles of those metal-line transitions from the system at z = 2.462 
without VPFIT solutions. The dashed vertical line corresponds to 2;=2.4628. 
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Fig. 4. — Velocity profiles and VPFIT solutions for the transitions from the system 2;=2.457. 
The dashed vertical line corresponds to Z—2A5785. Those features drawn in dotted marks 
are due to absorption from other systems. 
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Fig. 5. — Velocity profiles for the metal-line transitions for the system at 2;= 1.955. The 
dashed vertical line corresponds to 2=1. 9555. 
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Fig. 6. — Comparison of the low-ion profile Fe II 1608 with the Lya profile for the system at 
2;=!. 955. The overplotted curve in the Lya plot is a Voigt profile corresponding to 2;=1.9555 
and A^(HI) = 1.5 x lO^^cm-^. 
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Fig. 7. — Velocity profiles and VPFIT solutions from the system at z=2.325. The dashed 
vertical line corresponds to 2;=2.32612. 
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Fig. 8. — Velocity plot of the C IV 1550, Si II 1527 and Lja transitions for the system at 
z— 2.325. The transitions all span nearly the same velocity interval and appear consistent 
with profiles of Lyman limit systems. 
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Fig. 9. — Velocity profiles of two Mg II metal systems at (a) Z=1A1Q and (b) 2;=!. 699. The 
dashed vertical lines are at 2;=!. 47607 and 2=1.699 respectively. 
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Fig. 10. — Apparent column density [iVa(-;;)] for Fe II in the (a) ^=1.476 and (b) 
2;= 1.699 systems. In (a), the regions where Fe II 2374 dominates indicate hidden saturated 
components. On the other hand, the features in (b) where the Na{v) profiles do not coincide 
are due to blends with other absorption systems. 
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Fig. 11. — Apparent column density [iV(j(u)] profiles for Ni II in the z=2A62 system. The 
good agreement between the two profiles suggests there are no hidden saturated components. 
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Fig. 12. — Apparent column density [iV<j(v)] profiles for Cr II in the 2;=2.462 system. The 
Cr II 2062 profile is clearly stronger over the entire velocity interval. We suggest this is due 
to blending in the Cr II 2062 profile by Zn II 2062. In particular, we believe the features at 
—60, 40 and 60 km s~^ are due solely to Zn II absorption. 
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Fig. 13. — Velocity plot of the Fe II 1608 and Lya transitions for the z=2A57 and 
^=2.462 systems. The Voigt profile overplotted on the Lya velocity plot was derived with 
23 individual HI components (See Table 10). 
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Fig. 14. — (a) Fractional ionization versus log[A^(HI)] of a constant density (n = 0.1 cm~^), 
plane-parallel Hydrogen disk assuming a temperature of 10'^ K. The ionizing spectra is an 
attenuated power law spectrum, calculated by Madau (1992) at an average redshift of z=2A6. 
(b) Fractional ionization versus disk depth (expressed in cell number where cell 1 and 100 
are the outer faces) for a Hydrogen system with log[iV(HI)] = 20.4. Note only the edgemost 
cells are significantly ionized. 
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Fig. 15.— 



Comparison of our ionization model with CLOUDY simulations. 
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Fig. 16. — CLOUDY simulations for a series of systems with varying hydrogen volume 
density assuming [Z/H] = —0.5 and requiring log[iV(HI)] = 20.38. The ionizing spectra is 
an attenuated power law spectrum, calculated by Madau (1992) at an average redshift of 
^=2.46. 
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Fig. 17. — Plot of the relative logarithmic column densities of (a) Si+ vs. Si^"*" and (b) Al"*" 
vs. A1++. The [A1"'"/A1"^+] ratio was derived assuming logA^(Al+) = 12.0 at each pixel and 
should be considered a lower limit. The two figures combined with our CLOUDY simulations 
further suggests that this system is predominantly neutral. 
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Fig. 18. — Relative cosmic abundances of Zn, Fe, Ni, Si, and Cr versus Zn for 3 velocity 
features of the Z—2A62 system. Note how the relative abundances are essentially the same 
over all 3 features. 
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Fig. 19. — Logarithmic abundances of the low-ions Fe, Cr, Si, Ni, Zn relative to standard 
solar abundances versus condensation temperatures for the z=2A62 system (solid dots) and 
for the line of sight in the ISM toward ( Oph (open squares). Note that the Zn abundance 
for the damped Lya system is a lower limit whose value is most likely ~ —0.26. 
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Fig. 20. — Percentage of total column densities in 5 features of the z=2.462 system for the 
Fe II, Ni II, Si II, Cr II, and Si IV transitions. Note the obvious differences between the 
low and high-ions. The errors were derived from the apparent optical depth method and 
therefore should be taken only as an approximation. 



